Quantum Phase Transition in a Multi— Level Dot 
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We discuss electronic transport through a lateral quantum dot close to the singlet-triplet degener- 
acy in the case of a single conduction channel per lead. By applying the Numerical Renormalization 
Group, we obtain rigorous results for the linear conductance and the density of states. A new 
quantum phase transition of the Kosterlitz-Thouless type is found, with an exponentially small en- 
ergy scale T* close to the degeneracy point. Below T* , the conductance is strongly suppressed, 
corresponding to a universal dip in the density of states. This explains recent transport measure- 
ments. 
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Introduction. In the past years semiconductor quan- 
tum dots have gained considerable attention as tunable 
magnetic impurities jjj . Due to their small size electronic 
transport is strongly influenced by Coulomb blockade j| . 
A well-known many-body phenomenon, the Kondo ef- 
fect, was found in quantum dots j| with odd electron 
number, as predicted earlier jj]. In these systems, a sin- 
gle unpaired spin is screened at low temperatures, giving 
rise to an enhanced conductance at low bias. 

Remarkably, a similar conductance enhancement was 
recently also observed for an even number of electrons 
both in vertical || and lateral ||] GaAs quantum dots. 
This can be understood by taking into account the strong 
intra-dot electronic exchange coupling which is ferro- 
magnetic, similar to Hund's rule in atomic physics. Tun- 
ing of the level spacing by an external magnetic field then 
induces a singlet-triplet transition in the ground state. 
This additional degeneracy leads to the enhanced con- 
ductance at low temperature P-fTPf . 

Model calculations for the singlet-triplet transition, 
except [ |l0[ , have so far mainly considered the case where 
the orbital quantum number characterizing the levels in 
the dot is also present in the leads. In particular, two con- 
duction channels per lead have been taken into account. 
While this is appropriate for vertical devices ||, recent 
measurements on lateral quantum dots Jll]] suggest an in- 
terpretation in terms of a single conduction channel per 
lead, i.e. strong orbital mixing. In the following, we will 
focus on this situation. 

Our main tool of analysis is Wilson's Numerical Renor- 
malization Group (NRG) jL2j, a non-perturbative ap- 
proach to quantum impurity systems. In contrast to 
mean-field or scaling calculations, this method does not 
rely on any assumptions regarding the ground state or 
leading divergent couplings. We will find that this is cru- 
cial in the present analysis. 

The Model. We consider a two-level Anderson impu- 
rity model as shown in fig. |l|. The Hamiltonian H = Hr,+ 
H R + H D + H T contains two leads H r = J^kr e kr a \ ar a kar 
with r = L/R. The isolated dot is described by 



Hr 



J~] e dn dl-d„ 



JS 1 S 2 + E c (N-Af) 2 (1) 



where n = 1, 2 denotes the two levels, N = Y^na dn a d na 
is the total number of electrons occupying the dot and 
S„ = (1/2) Y,oo> dna^a<y'd nry , are the spins of the two 
levels. Furthermore, we have introduced the charging 
energy Ec and an exchange coupling J which arises due 
to Hund's rule. We choose J\f = 2 in order to achieve 
double occupancy of the dot. Through the energies e^n 
the level spacing Ae = edi — e<ii as well as the precise 
position in the Coulomb blockade valley can be tuned. 
Experimentally, M depends on the gate voltage, while 
Ae is controlled by an external magnetic field. As a 
consequence of Hund's rule, the intra-dot exchange is 
ferromagnetic (J < 0). Therefore, for Ae = —J/4 the 
three triplet configurations |1, 1 >= dj|G?2j 

|0 >, 1 1, > = 

(l/v^)(4 T 4j. +4|4tI° >' I 1 ' -1 >= d U d 2il° > and 
the singlet 1 0, >= d^dJjjO > are degenerate. Moti- 
vated by the small g-factor in GaAs [[|, we neglect the 
Zeeman splitting of the triplet states. 




FIG. 1. Two-level quantum dot (|l|) with single-channel 
leads (chemical potentials fir, and fJ-n). V denotes the tun- 
neling matrix element between the dot levels and the leads. 

Finally, tunneling between the dot and leads is mod- 
eled by H T = Y,kanr( V nral a d n<7 +h.c.) where we neglect 
the energy dependence of the tunneling matrix elements 
V nr and in the following take them to be symmetric and 
identical for both levels, that is V nr = V. The intrinsic 
line width of the dot levels due to tunnel coupling to the 



1 



leads is T = T L + T R with Y L/R = 2ir\V\ 2 N L/R where 
Nl ir is the density of states in the leads. This model has 
been studied before (l3) without discovering the quantum 
phase transition we describe in this paper. 

Transmission probability. We are interested in calcu- 
lating electronic transport through the dot ([!]) close to 
the singlet-triplet transition. To this end, we use the 
generalized Landauer formula Il4| 



^ J duj (f(uj - hl) 



f{u) - fi R )) T(uj) 



(2) 



with the Fermi function f(u) and the transmission coef- 
ficient 



TV) = - E 



L~nR 



n,n' ,<r 



T L +T R 



ImG m - a (w). 



(3) 



Here we have introduced the retarded dot Green's func- 
tions G nn > a {t) = -i9{t) < {d n<T (t),dl, a } >. In the fol- 
lowing we focus on the low bias regime, where T(u)) can 
be evaluated in equilibrium, using the Numerical Renor- 
malization Group. For a detailed description of this tech- 
nique see Ref. [l2[| . Note that the equilibrium transmis- 
sion T(w) also yields an approximation to the differential 
conductance dl/dV at finite bias. 



Ae/D = 0.3. Note the unusual shape of this peak, which 
is due to the "underscreening" of the local moment - a 
free spin 1/2 remains present in the ground state |l~5| ] 
and leads to logarithmic corrections to Fermi-liquid be- 
haviour. Nevertheless, as for the spin 1/2 Kondo effect, 
the transmission reaches the unitary limit at low temper- 
atures. Due to systematic numerical errors in the NRG 
calculation, this limit is underestimated by about 10%. 

In the regime Ae/D > 0.5, both electrons occupy the 
lower dot level and the ground state is a singlet. Re- 
markably, in this case a pronounced dip arises within the 
Kondo peak, leading to a strongly reduced transmission 
at low energy. The residual value T(0) is independent 
of Ae and is determined by the position in the Coulomb 
blockade valley. In particular, it vanishes in the center of 
the valley where tdi = —e<i2- 

In order to demonstrate clearly the importance of the 
number of conductance channels, we show results for the 
case of two transmission channels onto the dot ([!]) in 
the inset of fig. ||. In this case, the local spin is always 
completely screened by the leads. One obtains a con- 
ventional Zeeman type splitting of the conductance peak 
due to the energy difference between singlet and triplet. 
The sharp dip described above does not occur here and 
is thus characteristic for the single channel situation. 




FIG. 2. Transmission coefficient at zero temperature for 
different level spacings Ae = e^i — £d2 at N = 2, Ec/D — 1, 
edi/D = 1, T = 0.57D and J/D = -2. The bandwidth is 
given by 2D. Inset: For comparison, we show the transmis- 
sion in the case of two conduction channels per lead, identical 
dot parameters, V = 0.2873 and Ae/D = 0.4, 0.5, 0.55, 0.6, 0.7 
(from top to bottom). 

In Fig. U we show results for the transmission as a func- 
tion of the level spacing Ae, corresponding to a variation 
of the external magnetic field in the experiment. Our 
unit is the half bandwidth D of the conduction electrons. 
For Ae/D < 0.5, both orbitals are equally occupied and 
the dot is in a triplet state (S — 1). Due to the hybridiza- 
tion with the leads, this local spin is partially screened, 
giving rise to the Kondo resonance shown in fig. ^| for 




FIG. 3. Scaling of the dip at zero temperature. Param- 
eters are chosen as in fig. El but with a smaller broadening 
r = 0.25-D. Note the "pinning" of the transmission at uj = 0. 
Inset: Rescaled dip T(lo/ui*) for parameter values close to 
the transition. The curves collapse onto a universal scaling 
function with a single parameter uj* (dip width). 

We find that for small level spacing the dip has a uni- 
versal scaling form which is completely determined by its 
width uj* . This can be seen more clearly in fig. |3[ where 
we focus on the singlet side close to the transition. The 
scaling curve is extracted in the inset. 

Linear conductance. Using the current formula (j^), 
we now determine the behaviour of the conductance as 
a function of temperature. Results are shown in fig. ^. 
In the triplet regime, upon lowering of the temperature, 
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the conductance rises monotonously up to the unitary 
limit due to the partial screening of the local spin S — 
1. Note that the associated Kondo temperature Tk is 
extremely low. This may be the reason why the triplet 
Kondo effect has so far not been observed experimentally, 
though some indications have been seen in Ref. f§. Close 
to the singlet-triplet transition, Tr is strongly enhanced. 
On the singlet side, we find a "bump" type behaviour of 
the conductance when T is lowered: After an initial rise 
due to the Kondo effect, G(T) decreases strongly at T < 
T* . with a small residual value for T — > determined by 
the position in the Coulomb blockade valley. Note that, 
like the increase, the decrease of G(T) is logarithmic, 
indicating a two-stage Kondo effect. In particular, the 
T — > behaviour of G(T) is again universal and can be 
characterized by a single fit parameter T* ~ to* . 





<x> Ae/D = 0.4 
□-□ 0.46 
o ~o 0.487 
a-a 0.49 
*~* 0.495 




0.0001 



0.01 



T/D 



FIG. 4. Linear conductance for identical parameters as 
in fig. ^ in three different regimes: on the triplet side 
(Ae/D = 0.4,0.46), at the transition (Ae/D = 0.487) and 
on the singlet side (Ae/D = 0.49, 0.495). 

Quantum phase transition. Here we present a physical 
explanation of the above results. It has been suggested 
earlier |l(| that the singlet-triplet degeneracy of the dot 
can be parametrized in terms of a two-spin S = 1/2 
Kondo model. Formally, this is achieved by performing 
a Schrieffer-Wolff projection [|L7j of our two-level An- 
derson Hamiltonian on the (almost) degenerate subspace 
spanned by the four states |0, >, |1,1 >, |1,0 > and 
1 1,-1 >• One obtains an effective Hamiltonian of the 
following form 



H eS = H L +H R + J x S lS + J 2 S 2 s + /SiS 2 . 



(4) 



where both Kondo spins are coupled to the same conduc- 
tion channel s = (l/2iV)E 



A- A- 



/ a; 



i aw. 



Here, N is 



the number of fc-states in the leads and we have already 
taken a symmetric combination of left and right lead ac- 
cording to Ofccr = {akaL+a-kaR) / v2. Additional potential 
scattering terms have been neglected. Note that the Si, 
S2 introduced above are fictitious spins, different from 
the original levels. 



To leading order, the parameters in (Q) are given by 
the following expressions (note that V ~ 1/yN): 



J m = 2NV' 



1TV2 



1 



f-di + E C - J/4 e d2 + E C - J/4 
1 1TV2 \ 



(5) 



e d i -E c + J/4 e d2 -E c + Jj '4/ 



I = 2NV 1 



1 



f 



v e d i - E c + J/4 e d2 - E c + J/4 
J 



£d2 - Ec 



(6) 



In particular, we find J\ 7^ J 2 . The effective direct ex- 
change I will be a function of the level splitting Ae in 
the original model. 

The Hamiltonian (|4|) has been analyzed recently Jl8| . 
Depending on the strength of /, the ground state of the 
two spins is either an inter-impurity-singlet or a triplet. 
The associated transition at I = / cr it is of the Kosterlitz- 
Thouless type. The triplet side corresponds to an under- 
screened S = 1 Kondo model, while on the singlet side, 
a two-stage screening process of the two impurities has 
been found for small AI = I — I cr n > 0. First, the larger 
one of the two couplings (e.g. Ji) leads to a screening 
of the corresponding spin Si by the Kondo effect, thus 
decoupling S2 from the conduction band for T < Tk. 
The effective inter-impurity exchange AI then leads to 
a second Kondo effect for S>2. At low temperatures, the 
two spins form a singlet with a binding energy 



T* ~ exp(-TV/AI) 



(7) 



that is indeed exponentially small in the distance A7 « 
Ae — Ae cr it from the critical point. This argument holds 
only as long as A J < Tk- For larger values of the effective 
exchange, AI provides a cutoff on the "first" Kondo effect 
and the singlet binding energy is then linear, T* « AJ, 
as a function of the level splitting. In both cases, trans- 
port at T < T* is strongly suppressed due to the singlet 
formation which leads to the dip in the density of states. 

In order to demonstrate that this is the correct low- 
temperature scenario of the two-level quantum dot, we 
have calculated the characteristic energy scale T* deter- 
mined by the width of the dip. This calculation has been 
performed for the full dot Hamiltonian (|l|). In fig. ||, 
we give T* as a function of the distance from the criti- 
cal point. Clearly, for large |Ae — Ae cr ;t| the dip scales 
linearly, while close to the critical point a crossover to 
an exponential dependence occurs. Note that different 
level broadenings lead to largely different Kondo tem- 
peratures; in particular, for T/D = 0.063 only the linear 
behaviour of T* is seen because Tk is extremely small. 

At this point we would like to point out the robust- 
ness of our results with respect to parameter changes in 
the model. We have chosen our quantum dot to be at 
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a generic position in the Coulomb blockade valley, thus 
demonstrating that the quantum phase transition and 
the suppression of low-temperature transport discussed 
here are not restricted to special situations like particle- 
hole symmetry. The dip is also found when the broaden- 
ing of the two levels is tuned to different values and/or 
when an asymmetry between the coupling to the right 
and left lead is introduced. 
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FIG. 5. Scaling of the characteristic temperature scale T* 
(width of the dip) versus the distance from the transition 
point for E c /D = 1, e dl /D = 1 and J/D = -2. 

Conclusion and experimental relevance. Motivated by 
recent transport experiments we have studied trans- 
port through a lateral quantum dot modelled by two lev- 
els close to the Fermi surface coupled to a single conduc- 
tion channel in the leads. Correlations between different 
electrons in the dot are taken into account via the charg- 
ing energy and a ferromagnetic exchange coupling due to 
Hund's rule. 

At the associated singlet-triplet degeneracy, we find 
a quantum phase transition of the Kostcrlitz-Thouless 
type, as in the two-impurity model |l8[ ]. On the triplet 
side of the transition the conductance simply increases up 
to the unitary limit upon lowering of T. On the singlet 
side, we find a non-monotonic behaviour of the conduc- 
tance as a function of temperature ( "bump" ) correspond- 
ing to a characteristic "dip" in the transmission, which is 
also expected to be seen in the differential conductance. 
In particular, the width of the dip represents a new low- 
energy scale - the singlet binding energy - which becomes 
exponentially small close to the transition. For two con- 
duction channels, none of the two effects is observed. 

These findings are in good agreement with recent con- 
ductance measurements for a lateral quantum dot fi~l| l 
close to the singlet-triplet transition. In this system, or- 
bital symmetry is not conserved during tunneling. Both 
the non-monotonic behaviour of G(T) and the sharp dip 
in dl/dV have been observed, in contrast to previous 
studies of vertical quantum dots p[ . 



We therefore conclude that the number of conduction 
channels plays a crucial role for low-energy transport 
properties of a quantum dot. Symmetry and physical 
behaviour of such a device are thus strongly related. 
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